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Abstract We used a two-dimensional finite element 
method to analyze the load-slip characteristics of nailed 
wooden joints sheathed with a panel. We used tests of nail 
bending, nail shank embedment in a wood or a panel, nail- 
head embedment in a face of a panel, nail withdrawing from 
a wood, friction between a wood and a panel, and initial 
axial forces of the nails. The values of the material proper¬ 
ties for the analytical models were determined with the 
above tests. With a conventional one-dimensional analysis 
it is impossible to evaluate the shearing performance accu¬ 
rately because the axial forces of the nail are not calculated. 
Therefore, we used two-dimensional geometric nonlinear 
analysis. The computed load-slip curves closely matched the 
original experimental results; and when the friction and 
initial axial forces were considered, the computed curves 
were more likely to match the original experimental results. 
Our results suggested that it is better to use geometric non¬ 
linear analysis to analyze the shearing characteristics of 
nailed joints, and that friction and axial forces are important 
factors for accurate analysis. 
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Introduction 

Analytical methods to evaluate load-slip characteristics of 
mechanical wood-joints have been established in various 
studies. - Generally, the yield load is estimated by a 
method based on the European yield theory, the stiffness is 
estimated by a method based on the theory of beams on 
elastic foundations, and the load-slip curve is calculated by 
methods based on the theory. 

We believe there are two problems with these methods. 
One is that the calculated load is smaller than the actual 
experimental load owing to not considering the lack of ef¬ 
fect of the fastener’s axial forces, which are the forces that 
restrain axial displacement of the fasteners by a nut or a 
nail-head. It is impossible to calculate the fasteners’ axial 
forces by conventional one-dimensional methods. The 
other problem is that the calculated stiffness is also quite 
different from the actual experimental one. 10 

The purpose of this study was to find a better way to 
calculate the load-slip curve of the nailed wood joints 
sheathed with a panel. We tried various methods, including 
one based on the small deformation theory and geometric 
nonlinear analysis in models of the nailed wood joints, tak¬ 
ing into consideration the friction force between the wood 
and the panel and the initial axial forces of the nail. 

An analytical method that considers the axial and fric¬ 
tion forces was presented by Hirai and Wakashima. 68 We 
think that with the geometric nonlinear analysis we used the 
axial forces can be estimated more accurately and directly 
by calculation, and the fastener’s large deformation can be 
predicted more exactly by the convergence calculation than 
when using Hirai’s method. 


Materials and methods 

Materials 

We used medium-density fiberboard (MDF) (thickness av¬ 
erage 8.95 mm, density average 0.73 g/cm 3 ) and conifer lum- 
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Fig. 1. Single shearing tests of 
nailed wood joints. Nail heads of 
type 1 touched the face of the 
medium-density hberboard; 
(MDF). Type 2 nails were 
pounded more deeply than type 1 
nails based on the thickness of 
the nail head 
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Fig. 3. Nail pull-out tests. Displacement of the crosshead was mea¬ 
sured. S-P-F , conifer lumber for wood frame construction 


ber for wood frame construction (S-P-F 204) (density aver¬ 
age 0.54 g/cm 3 ). The nails were CN50. 


Methods 

We constructed the test specimens for single shearing tests 
of nailed wood joints, among which there were two types 
(for type 1 the nail head was not pounded into the panel, 
and for type 2 the nail head was pounded into the panel) 
(Fig. 1): test specimens for nail-bending tests (Fig. 2), test 
specimens for nail pull-out tests (Fig. 3); test specimens for 
lateral nail resistance tests (Fig. 4); and test specimens for 
nail head pull-through tests (Fig. 5). For the lateral nail 
resistance tests, the gaps between the specimen and the jig 
were within 0.02 cm, and the length/diameter ratio was 3.21. 
These tests, except for the nail head pull-through tests were 
done under monotonic loading. 

Figure 6 shows a typical nail head pull-through test. First, 
the nail is not fully pounded into the MDF, and the nail 
head does not touch the face of the MDF. The nail is then 
pulled at a steady speed of 2mm/min. As soon as the nail 
head touches the face of the MDF, the pulling is stopped 
and the force is fully released. Pulling begins again at the 
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Fig. 4. Lateral nail resistance tests 
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Cross-head speed 
: 2mm/min 




X-X’ cross-section 


MDF 

(9*75*75 mm) 
n~7 


Fig. 5. Nail head pull-through tests 




Loading Removal Loading 
of load 


Fig. 6. Loading during the nail head pull-through test 


An element after load application 



An element before load application 

Fig. 8. Difference between the small deformation theory and geomet¬ 
ric nonlinear analysis. Left Small displacement theory. Right Geomet¬ 
ric nonlinear analysis. P xl , P yl , M u forces at the first node; P x2 , P y2 , M 2 , 
forces at the second node; a, slope of the element 


cell touching the S-P-F. The coefficient was calculated by 
the formula 

c df = P f /(w dl + w spf ) (i) 

where C df is the coefficient of dynamic friction; P f is the 
friction force as indicated in Fig. 7; W dl is the weight of the 
dead load (N); and W spf is the weight of the S-P-F (N). 


Numerical solution 

We simulated the load-slip characteristics of the single 
shear nailed wood joints sheathed with MDF. The models 
were computed by the load-incremental method of the fi¬ 
nite element method (FEM). In addition, the stiffness and 
strength were calculated by conventional analytical meth¬ 
ods based on the theory of beams on elastic foundations and 
the European yield theory. 



Estimated speed: 2mm/min 

Fig. 7. Rubbing test and valuation method for coefficient of dynamic 
friction. P f , dynamic friction force 

same speed until the force decreases to 80% of the maxi¬ 
mum load. 

Figure 7 shows the schema for the rubbing test and for 
evaluating the dynamic friction coefficient. The MDF was 
moved in a horizontal direction, and the friction forces be¬ 
tween the MDF and the S-P-F were measured with the load 


Finite element method 
Techniques 

The FEM is classified into two categories. One is a method 
based on the small deformation theory, and the other is 
geometric nonlinear analysis for large deformation. With 
the small deformation theory, the equilibrium equations are 
based on the undeformed geometry before load application. 
For the geometric nonlinear analysis, equilibrium is based 
on the real deformed configurations, as shown in Fig. 8. The 
details are described elsewhere. 11,12 

When the load on a geometrically nonlinear structure is 
applied in small increments, the force and displacement 
may be related by 

i s U AD L = i AF ) m ( 2 ) 

where {AF} m and {AD} m are increments of forces and dis¬ 
placements. The tangent stiffness matrix [S t ] m changes with 
the change in geometry, and [S t ] m depends on the internal 
forces. [S t ] m is expressed as the sum of two matrices 

LL = NL„ + 


( 3 ) 
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where [S e ] m is the elastic stiffness matrix, and [S g ] w is the 
geometric stiffness matrix. The two matrices vary during the 
analysis. If the change in geometry and internal forces are 
neglected, [S t ] m is the tangent matrix based on the small 
deformation theory. 

In our geometric nonlinear analysis, we considered only 
the geometric nonlinearity of the beam elements and in¬ 
tentionally did not consider that of the axial springs. The 
material nonlinearity was considered in both the small de¬ 
formation theory and the geometric nonlinear analyses. As 
for the beam elements in the geometric nonlinear analysis, 
[S.L was expressed as a function of length and slope a of 
the beam element internal forces P% in Fig. 9, and Young’s 
modulus of the fastener. [S g ] m was expressed as a function of 
length and slope a of the beam element and internal forces 
P% and P *i i n Fig- 9. As for the beam elements in the small 
deformation theory, the length and slope of the beam ele¬ 
ment were assumed to be constant in the initial position, 
and the internal forces were neglected. Then [S g ] m was zero 


Pyi 



Fig. 9. Coordinate system for a beam element. Px%, Py *, Mf, internal 
forces at the first node in local directions; Px f, Py f, Aff, internal 
forces at the second node in local directions; a, slope of the beam 
element 


matrix, and [S e ] m was expressed as a function of Young’s 
modulus of the fastener. As for the spring elements in both 
analytical methods: [S t ] m was based on the small deforma¬ 
tion theory, [S g ] m was zero matrix, and [S e ] m was expressed 
as a function of the spring constant. The stiffness matrices 
were described previously. 13 

Analysis involving iteration is necessary if both material 
and geometric nonlinearities are considered. The analysis 
was done using Newton-Raphson’s technique. 

The iteration cycles were terminated when the following 
criterion was satisfied 

(w-wr=< 8 ({ f } 7 { f }f (4) 

where fi was the tolerance value, which we made 0.005 for 
the small deformation theory or 0.01 for the geometric non¬ 
linear analysis because the process of convergence in the 
geometric nonlinear analysis was more complicated than 
that in the small deformation theory. This criterion ensured 
that the out-of-balance forces {g} were small compared with 
the incremental forces {F}. 

The matrix was calculated by the Gauss-Jordan elimina¬ 
tion, which included complete pivoting. 

Analytical models 

Table 1 shows scheme of the computation models. We made 
two beam-spring models for single shearing analysis of the 
nailed wood joints sheathed with MDF: model 1 and model 
2 (Fig. 10). The springs used are indicated in Table 2. 

In both models 1 and 2 the nail shank is modeled using 50 
beam elements with rotation spring elements at both ends 
of each beam element (Fig. 11), among which 20 on the side 


Table 1. Scheme of the computation models 


Model 

name 

Model 

no. 

Geometric nonlinearity 
(only the beam elements) 

Friction and 
initial axial force 

Object (type) 

Hl-sma 

1 

Neglected 

Neglected 

1,2 

Hl-geo 

1 

Considered 

Neglected 

1,2 

Sl-sma 

2 

Neglected 

Neglected 

1,2 

Sl-geo 

2 

Considered 

Neglected 

1,2 

Fl-geo 

2 

Considered 

Considered 

1 

F2-geo 

2 

Considered 

Considered 

2 


Table 2. Springs 

in the computation models 


Spring no. a 

Type 

Direction 

Experiment and object 

XI 

Axial 

X 

Pull-out tests of nails 

Withdrawing of nail from the S-P-F 

X2-X42 

Axial 

X 

Nail-head pull-through tests 

Embedment of nail-head in the MDF 

Y1-Y31 

Axial 

Y 

Lateral nail resistance tests of the S-P-F 
Embedment of nail shank in the S-P-F 

Y32-Y52 

Axial 

Y 

Lateral nail resistance tests of the MDF 
Embedment of nail shank in the MDF 

R1-R100 

Rotation 

6 

Bending tests of nail shank 

Nonlinearity of nail bending 


S-P-F, conifer lumber for wood frame construction; MDF, medium-density hberboard 
a Spring numbers are indicated in Fig. 9 
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Axial spring 
(Y-direction) 



Fig. 10. Model of the nailed joint sheathed with MDF. Squares, nodes; Table 2. We considered only the geometric nonlinearity of the beam 

broken lines, axial spring elements; heavy solid lines, beam elements. elements in the geometric nonlinear analysis 

XI to X42, Y1 to Y52, and R1 to R100 are the spring numbers listed in 


Rotation spring elements 



Fig. 11. Elements of the nail shank. E, Young’s modulus; I, geometri¬ 
cal moment of inertia; A, Cross-sectional area; L, length of the beam 
element; R1 and R2, nonlinear rotation springs (moment-rotation 
angle curves of R1 and R2 were assumed to be the same) 


of the nail head are in the MDF and the other 30 are in the 
S-P-F. The length of the beam elements of the models for 
type 1 is different from that for type 2. In the analysis F2- 
geo of type 2 (Table 1), the length of the part of the nail 
shank touching the MDF is 7.75 mm (8.95 — 1.20 mm) and 
that touching the S-P-F is 42.25mm (50.00 - 7.75mm), 
where 8.95 mm is the thickness of the MDF and 1.20 mm is 
the thickness of the nail head. The other analysis, the length 
of the part of the nail shank touching the MDF, is 8.95 mm 
and that touching the S-P-F is 41.05 mm. Nonlinearity asso¬ 
ciated with the axial forces of the nail is intentionally ne¬ 
glected, so each beam element can be assumed to be an 
elastic body. Nonlinearity associated with bending of the 
nail is simulated by the rotation springs (Fig. 12). The em¬ 
bedment is expressed by axial springs (Fig. 10). The forces 
of the axial springs were assumed to depend on the relative 
displacement in the x-direction or y-direction between the 
first and second nodes as described in Table 2. The broken 



9 0 0 


Fig. 12. Modeling method of nonlinearity of the nail shank. I, bending 
of the nail; II, beam element; III, rotation springs 

lines representing axial springs in Fig. 10 indicate the loca¬ 
tion of the first and second nodes of the spring elements; 
they do not express vectors of the springs. 

With conventional methods, the nail shank embedment 
is modeled by elements based on the theory of beams on 
elastic foundations. Theoretically, the calculated load-slip 
relation of the models using fewer axial springs is different 
from that of models using elastic foundations, but the differ¬ 
ence can be reduced to a negligible level by increasing the 
number of axial springs. 

In model 1 the nail head is modeled using 40 beam ele¬ 
ments, each of which is assumed to be an elastic body. In 
model 2 the nail head is not modeled and is assumed to have 
free rotation. 

Input data processing 

The original experimental data were processed and mani¬ 
fested by polygonal lines before the FEM analysis (Fig. 13). 
The load-displacement relation of the axial springs and 
Young’s modulus of the beam elements were determined 
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Fig. 13. Input data of the springs for the computation models 


i 



CL 20mm 


Fig. 14. Beam-spring model for nail-bending tests. P/2, half of the 
applied load in the nail bending test; CL, center line 

from the average of the experimental values. Young’s 
modulus for the beam elements is 185.5 GPa. 

The moment-rotation relation of the rotation springs is 
determined by trial and error so that the calculated result of 
the beam-spring model (Fig. 14) fits the load-displacement 
curve of the nail-bending tests (Fig. 2). Only the right half of 
a nail-bending test is simulated in the model shown in Fig. 
14 because the whole test body is symmetrical. 

It is assumed that the axial springs from X2 to X48 in Fig. 
10 show resistance against the deformation only in the area 
where the nail head is embedded in the MDF. For the 
computation models HI and SI in which the friction be¬ 
tween MDF and S-P-F is not considered, the curve segment 
from 2 to 5 in Fig. 6 was used as the input data. For Fl-geo 
and F2-geo, we considered not only the friction but also the 
initial forces. For Fl-geo, the curve segment from 3 to 5 was 
used as the input data; and the tensile force, which was 
equal to P3 was applied to the nail head and the nail tip 
because an initial force P3 was generated when the nail 
head touched the MDF. 

For F2-geo, the curve segment from 4 to 5 was used as 
the input data; and a tensile force, which is equal to the 
initial force PA, was applied to the nail because the nail 
heads of type 2 were pounded into the MDF (the nail head 
thickness was 0.12 cm, and the nail head surface was at the 
same level as the MDF surface). 

The dynamic friction coefficient between the MDF and 
the S-P-F was 0.166, which was a calculated average value 
from the rubbing tests shown in Fig. 7. The friction force 
was calculated from the axial force of the spring X2 in Fig. 
10 . 


Conventional analytical methods 

The initial stiffness was calculated based on the theory of 
beams on elastic foundations. The yield load and the maxi¬ 
mum load were calculated by the European yield theory. 
The calculation was conducted referring to other people’s 
methods, and the parameters were similarly determined by 
lateral nail resistance tests and nail bending tests. 9 The 
maximum load was calculated from the maximum stress in 
the lateral nail resistance tests and the maximum moment in 
the nail-bending tests. The slope between 10% and 40% of 
the maximum load in the lateral nail resistance tests was 
used to calculate the initial stiffness. 


Results and discussion 

Figure 15 shows the results of the nail deflection calculated 
by Sl-sma and Fl-geo computational models. The right- 
hand side shows the results of the Fl-geo geometric nonlin¬ 
ear analysis. The figure shows that there were nail head 
embedment on the upper right side and nail withdrawal on 
the lower left side. In this case, the calculated nail length 
was near its actual length. The left-hand side shows the 
calculated results by the small deformation theory, in which 
the nail nodes did not move in the direction perpendicular 
to that of loading, there was no nail-head embedment or 
nail withdrawal, and the calculated nail length was longer 
than its actual length. This is because the axial forces were 
not counted using the small deformation theory in which 
the equilibrium is based on the node position before load 
application. 

Figure 16 compares the original experimental results and 
the calculated results, neglecting the friction force in the 
single shearing tests. Although the maximum loads of types 
1 and 2 were almost the same, their load-slip curves were 
different until the maximum load. The stiffness of type 1 
was less than that of type 2. 

The calculated curves by geometric nonlinear analysis 
were more consistent with the experimental curves than 
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Fig. 15. Nail deflected shape (theoretical) 
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Fig. 16. Comparison of the original experimental results of the single shearing tests and the calculated results, with the friction and initial force 
neglected 


those calculated according to the small deformation theory. 
The determining factors of the calculated ultimate shearing 
strength by each analytical method were different. The de¬ 
termining factors in the small deformation theory were 
based on nail bending, and those in the geometric nonlinear 
analysis were based on nail shank embedment in the MDF, 
as shown in Fig. 15. Therefore, the calculated curves by each 
method were different. 

With the small deformation theory the nail started to 
yield at loads below 0.8 kN, but with the geometric nonlin¬ 
ear analysis the yielding load exceeded 0.8 kN because of 
the effect of the axial force of the nail shank. The mecha¬ 
nism of this phenomenon is that the nail head resists em¬ 
bedment into the face of the MDF. Moreover, tensile stress 
is generated in the nail shank, which makes the nail shank 
difficult to bend and accounts for the load increase in the 
single shearing tests of the nailed joints. 

The results calculated by both the European yield theory 
and the theory of beams on elastic foundations are also 


shown in Fig. 16, in which K is the initial stiffness, P y is the 
yield load, and P max is the maximum load. The results calcu¬ 
lated according to the European yield theory were similar 
to those found by the small deformation theory but were 
different from the original experimental values. 

As for the initial stiffness, all the calculated values were 
almost the same but were different from the original experi¬ 
mental values. Moreover, the calculated stiffness by the 
model HI-geo in which the nail head was modeled was 
almost the same as that found by the model Sl-geo in which 
the nail head was not modeled. This indicates that the rota¬ 
tion of the nail head was hardly restrained. 

Figure 17 compares the experimental and calculated re¬ 
sults of type 1 with the friction and the initial axial forces 
considered. Figure 18 compares the experimental and calcu¬ 
lated results for type 2, with the friction and initial axial 
forces considered. 

For both types 1 and 2, the calculated load-slip curves, 
taking into consideration the friction and the axial force, 
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Slip (cm) 

Fig. 17. Comparison of experimental and calculated results, with the 
friction and the initial force considered: TYPE-1 



Slip (cm) 

Fig. 18. Comparison between experimental and calculated results, 
with the friction and the initial force considered: TYPE-2 


were overall similar to the original experimental ones. The 
friction and the axial force may reflect a difference of the 
nail head position or embedment depth between types 1 
and 2. 

As for the initial stiffness when considering friction, the 
calculated value closely matched the original experimental 
value, though not exactly. The difference between the cal¬ 
culated value and the original experimental value is due to 
the fact that in both of these analytical models the effect of 
nail bending in the lateral nail resistance tests was not ex¬ 
cluded and the friction coefficient was set up as a constant. 
We think that the difference will diminish if we get rid of the 
effect of the nail bending and adjust the friction coefficient 
depending on the displacement. 


Conclusions 

The results of our study suggest that geometric nonlinear 
analysis is a relatively better method for evaluating 
load-slip curves of nailed wood joints. The calculated 
results closely matched the original experimental ones 
when the friction between the S-P-F and the MDF as well as 
the initial axial forces of the nail are considered in the 
calculation. 

Furthermore, this analytical method may make it pos¬ 
sible to evaluate precisely the load-slip characteristics of 
various nailed wood joints by simple tests such as those we 
used in our study. They can be applied not only to the 
evaluation of structural designs but to new material devel¬ 
opment, such as with materials for structural panels. 

Acknowledgment The authors thank Noda Corp. for kindly offering 
us the MDF. 
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